Domain wall formation and spin reorientation in finite-size magnetic systems 
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We investigate the formation of stable one-dimensional Neel walls in a ferromagnetic slab with 
finite thickness and finite width. Taking into account the dipolar, the exchange and the uniaxial 
anisotropic crystalline field interactions, we derive an approximative analytical self-consistent ex- 
pression that gives the wall width in terms of ratios between the three different energy scales of the 
problem. We also show that, even when the crystalline anisotropy does not favour the formation of 
domain walls, they can yet be formed due to the dipolar interaction and the finiteness of the system. 
Moreover, using a Stoner-Wohlfarth approach, we study the magnetization reorientation inside the 
domains under the action of an external magnetic field and obtain the respective hysteresis loops, 
showing that their shapes change from squared to inclined as the width of the slab varies. Finally, 
we discuss possible applications of this model to describe qualitatively some recent experimental 
data on thin films of MnAs grown over GaAs substrates. 

PACS numbers: 75.70.Ak ; 75.60.Ch ; 76.60.Es 



I. INTRODUCTION 

Magnetism in the micro and nanoscales is a source of 
promising technological advances in a very broad range of 
interests, from spintronics and quantum computation [l[ 
to biophysics and pharmacology. Nowadays, with the im- 
provement of experimental methods for growth and char- 
acterization of magnetic thin films, this particular class 
of mesoscopic (quasiVtwo-dimensional systems has been 
largely investigated 0, 0] • Particularly, one of the prop- 
erties of such films that has been calling more attention is 
their magnetic domain structures and their dependence 
on temperature, film thickness and applied magnetic field 

0, H 0, S| • Many of these systems present a competi- 
tion between a short-range, strong interaction (exchange) 
and a long-ranged, weak one (dipolar), from which it is 
expected the emergence of spatially modulated configu- 
rations 0, [13]. However, rarely these theoretical models 
consider possible size effects due to the finiteness of the 
film width; on the contrary, usually they describe the 
films as infinite plates, which is fair description of most 
experimental systems (see, for example, (lll.[l^.[T3.[l3|). 

One class of magnetic materials that cannot be de- 
scribed as an infinite plate and that has been recently 
subjected to deep experimental analysis are the thin 
films of MnAs grown over GaAs substrates (MnAs:GaAs) 

01. In bulk, MnAs exhibits a simultaneous abrupt first- 
order magnetic/structural transition from a ferromag- 
netic, hexagonal phase (a phase) to a paramagnetic, or- 
thorhombic one (J3 phase) flf| . However, MnAs:GaAs 
films do not show this abrupt transition; conversely, a 
large region of coexistence between a and j3 phases arises 
from C to 50 °C, characterized by the formation of 
periodic (a + /3) stripes of constant width. The rela- 
tive witdths of the a and (3 phases varies with temper- 
ature while the total (a + 0) width remains constant 
E3, US, EH US]- Hence, to understand the domain 
structure inside each ferromagnetic stripe, and how it 



varies with temperature, it is important to consider its 
finite width, specially because in most experimental stud- 
ies it is of the same order of its thickness. 

In this work, we apply the general method of energy 
minimization already used in other contexts [2ll . I22L |23 | 
(such as infinite plates, nanoparticles) to study the for- 
mation of stable unidimensional domain walls in a slab 
with finite thickness and finite width. The main purpose 
of such common procedure is to obtain an expression for 
the total energy that includes the different interactions 
contributions. In our case, there are three terms: the ex- 
change term (which tends to unfavour sharp walls); the 
uniaxial anisotropic crystalline term (which can favour 
or not the formation of sharp walls, depending on the 
easy axis of magnetization) and the dipolar term. This 
last one is rather important in finite systems, as already 
pointed out by others [T3|; we show that, in the system 
considered here, it is fundamental to form stable walls, 
specially when the crystalline anisotropy does not favour 
them. From this energy expression, it is possible to dis- 
cuss the different solutions for the wall width depending 
on the three scales of energy involved. 

We also generalize the energy expression to include do- 
mains whose main magnetization axis is tilted by an angle 
(j> with respect to the normal direction to the film. With 
such expression, it is possible to study how the magne- 
tization is reorientated under the action of an external 
magnetic field in the direction of the easy axis. This 
is achieved by calculating theoretical hysteresis curves 
through a method similar to the one proposed by Stoner 
and Wohlfarth [24]. This procedure does not take into 
account nucleation or pinning effects, but only the ro- 
tation of the domains and may lead to values of mag- 
netic coercive fields that are not exactly the measured 
ones. However, as it is a microscopic method, and not 
a phenomenological one, the main properties predicted 
are expected to be followed by a variety of experimental 
systems at least qualitatively. 
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Figure 1: General picture of the coordinate system used to 
describe the ferromagnetic slab. 



Here is an outline of the article: in Section 2, we pro- 
pose a general expression describing the magnetization 
corresponding to N domains whose walls width is a in a 
slab of thickness D and width d and calculate the corre- 
sponding total energy, obtaining an approximative ana- 
lytical equation. Next, we minimize it with respect to a 
and discuss various possible solutions for the wall width 
depending on the axis of magnetization (if it is the easy 
or the hard one) and the relationship among the three 
distinct energy scales. In Section 3, we generalize the 
previous procedure to inclined domains and discuss the 
different shapes of the hysteresis loops that describe spin 
reorientation (i.e., if they are squared or inclined and the 
values of the coercive fields). Section 4 is devoted to 
discuss a possible application of the model developed to 
understand some properties of MnAs thin films (partic- 
ularly, recent experimental hysteresis curves). Section 5 
contains the conclusions and final remarks of the work. 



II. ENERGY MINIMIZATION AND WALL 
WIDTH 



anisotropic system, we will not define yet which of the 
axis is the easy one. To consider the role of the slab 
thickness on the formation of these domains, we bound 
the magnetization to be nonzero only inside it: 

M = [M x {x)x + M z (x)z] 9(D/2 + z)6(D/2 - z) (1) 

where 8{x) is the Heaviside step function. It can be 
shown that this formulation would be equivalent to con- 
sider sharp domain walls in the upper and lower borders 
of the slab such that their width is much smaller than the 
width of the walls along the x direction. As we are con- 
cerned with the latter, there is no significant role played 
by the former in what follows, and we can do a z — » 0. 

To investigate the formation of one-dimensional sharp 
domain walls along the x direction, we propose to study 
the configuration in which the z component of the mag- 
netization is given by: 



N f-11 i_1 
M z (x) = M J2 K ' i<^<- 



x-d(i - l)/N 



V2a 



-erfc 



x - di/N 



V2a 



(2) 



where M is the saturation magnetization, erfc(a;) is the 
complementary error function, N is the number of do- 
mains inside the slab and a corresponds to the wall width 
along the x axis, which will be varied to minimize the 
total energy. This model resembles the one used in ref- 
erence [23j| ; the choice of using the complementary error 
function is due to its analytical properties that will allow 
us to obtain simple expressions for the wall width. 

As we are considering a crystalline anisotropy such that 
My = 0, we can obtain M x using the fact that the module 
of the total magnetization is constant and equal to M . 
Hence, if we meet the condition: 



Through all this paper, we will use the coordinate sys- 
tem shown in figure [TJ the z axis corresponds to the slab 
thickness (which we shall call D), the x axis, to the slab 
width (which we shall call d) and the y axis, to the slab 
length (considered infinite for our purposes). The origin 
of the axes is located at the middle point of one of the 
slab's faces, in a way that the x axis points in the direc- 
tion of the other parallel face. We consider that there 
is a strong crystalline field that does not allow the mag- 
netization to point in the y direction (xz spin model). 
This is the case for many experimental systems and, spe- 
cially, for the MnAs thin films. But we will postpone the 
discussion about realistic applications of the model until 
Section 4. 

Accordingly, the only kind of unidimensional domain 
walls that may be formed in this system are Neel walls. 
Initially, we choose the magnetization to lie along the 
z direction, in such a way that the walls are along the 
x direction. Although we are considering an uniaxial 



cr<Cd, (3) 

then we can obtain a simple approximate expression for 
M x in terms of Gaussian functions: 



e i^— . (4) 

i=l 

A typical domain configuration described by Q and 
(0]) is shown in figure [21 where the z and x components 
of the magnetization are shown as a function of the slab's 
width. 

As we discussed in the previous section, several works 
[g|. fTol] have shown that, in systems with competing inter- 
actions, a spatially modulated configuration is expected. 
One may interpret these configurations as "spread" do- 
main walls with sinusoidal domains; however, in this ar- 
ticle, we will investigate the formation of sharp walls, for 
which condition ([3]) is expected to be satisfied. 
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Figure 2: Components (a) M z and (b) M x of the magnetiza- 
tion described by |[2|) and for a four domain configuration 
along the reduced slab's width x/d. 



With the aid of expressions <[2j) and Q , we can calcu- 
late the total energy of a specific domain configuration. 
However, it is necessary to make an assumption about 
the microscopic nature of the system if we want to in- 
clude not only the "macroscopic" terms concerning the 
dipolar and the crystalline anisotropy interactions, but 
also the "microscopic" exchange energy. For long wave- 
lengths, one expects that the particular lattice structure 
(i.e., if it is cubic, hexagonal, etc) will not substantially 
change the qualitative physical properties derived for an- 
other kind of lattice. Hence, to simplify the calculations, 
we follow [id] and choose the system lattice to be cu- 
bic, with lattice parameter a and gyromagnetic factor g. 
Therefore, it is straightforward to relate the saturation 
magnetization to the microscopic parameters: 



M = 



9Pb 



(5) 



where ps is the Bohr's magneton. Now, it is possible to 
obtain the three different energy terms. The exchange 
term is obtained from the "classical" Heisenberg nearest- 
neighbours Hamiltonian: 



H exc = -— y^m- t ? 



j 

2a 



d A r 



Wm x (f) + Vm z (r) 



where we moved to the continuum limit of the lattice (see, 
for instance, [HI) and defined fa = M/M . Using {TJ, 
(121) and (@]), we obtain that the exchange energy density 
is given by: 



E, 



V 



2J 



2N + tt(N - 1) 
8y/r?{d/a)(a/a) 



(6) 



To obtain the dipolar energy, as there are no free cur- 
rents, we can use equation (see, for instance, the Magne- 
tostatics chapter of |26|): 



F _ 1 



(f)p(f)d 3 r, 



(7) 



where p = — V-M is the effective magnetic charge density 
and 4>(f) is the scalar magnetic potential, which satisfies 
the Poisson equation: 

V 2 = -At: p. 

Taking the Fourier transforms of p(f) and 4>{r), the 
Poisson equation can be easily solved and equation J7]) 
can be written as: 



Edip — 2tc 



p(k) 



k 2 



-d A k. 



Substituting expressions (!]), J2j) and a straight- 
forward calculation yields, for the magnetostatic energy 
density: 



E 



dip 



V 



1 1 
9 t*B 



+ e z ( -, N, p 



where p = d/D is the slab's aspect ratio and: 



(«) 



e '(l> N > 



2irp I du- 
lo 



- 2 (f) 2 



P 



sin 2 w/2 /u(iV-l)\ sinu/2 

1 2 cos ' v M 1 



sin^ u/2N 



2N 



sin u/2N _ 



-u/2p 



-l/< / due-^) 2 e -^sinh tan 2 (JL) [l - (-l)*cosu] . 



(9) 



Finally, we can calculate the uniaxial crystalline 
anisotropic term (see, for instance, the Magnetic 
Anisotropy chapter of [231]): 



E, 



cryst 



= -AK / tfVm 2 , 



4 



where AK, the anisotropy constant, can be positive or 
negative, depending if the x axis is the easy (AK > 0) or 
the hard one (AK < 0) . Evaluating this calculation, we 
obtain that the density of anisotropic crystalline energy 



E, 



cryst 



(N - 1)Jtt [o\ , x 

^ y^— y . do) 



Hence, the total energy is given by: 



Ek, 



E K 



E d 



If) 



E, 



cryst 



(ii) 



As we are assuming condition to be satisfied, we 
can make a further approximation to obtain a simpler 
expression for the total energy. Using such condition, we 
can approximate the exponential e~ u in the inte- 

grals ((9]) to 1, as long as we take an appropriate upper 
limit to them: 



(0,N,p) 



(l> N >> 



(12) 

This procedure implies in errors of the order of 10% to 
20%, if compared to numerical calculations. The choice 
of the upper limit of the integrals has a small impact 
on the final result since we are going to apply a self- 
consistent method in the end of the calculation. There- 
fore, the qualitative physical properties are still valid in 
this approximation, and we can minimize l[TTj) as: 



dE tl 



d 2 E t , 
da 2 



> 0. 



To solve these equations, it is convenient to use the 
following auxiliary variables: 



J 



a-' \9 2 V B l a 
(AKa 3 ) 3 



-1/3 



1/3 



/ 2N + n(N - 1) 
Uv5F(a/<0e* (0,N,p) 

2(N- 1) 3 tt 2 



-1/3 



1/3 



27 (ae x (0, N,p) /d) 2 (2N + n(N - 1)) 



(13) 



Then, the equation for the domain wall width can be 
written as: 



(14) 



a I \ a 



It is easy to see that the auxiliary variables l|13p are 
just relations between the three different energy scales 
involved in the system: in a', there is the ratio between 
the typical value of the exchange energy and the typi- 
cal value of the dipolar energy. The other term is just 
a numerical one and depends only on the ratio a/d, the 
number of domains N and the aspect ratio p. For thick- 
ness of the order of hundreds of lattice parameters, this 
numerical factor is usually of the order of 10 _1 , what im- 
plies that the relationship between the typical exchange 
and dipolar energies will determine the order of the wall 
width. From this, it is clear that when J ^> g 2 ^ijj/a 3 no 
sharp walls would be formed, as it would be expected. 

The parameter A is a relationship between the three 
types of energy and can be positive or negative, depend- 
ing if the x axis is the easy or the hard one, respectively. 
The numerical factor again depends only on a/d, N and 
p, and for thickness of the order of hundreds of lattice 
parameters, it is usually of the order of 10 _1 . 



Let us study the solutions of (JT4J) ; independently of 
the sign of A, equation lfl4|) always has only one positive 
solution. If A > 0, this solution is: 



— I = A 

a 



X 2 



2/3 



+ (l + Vl + A 3 ) , 
while, for — 1 < A < 0, we have: 



A + 



(1-Vl + A 5 ) 



2/3 



/ / \ 2 / 3 

+ (l - v/T+a^) , 

and in the case where A < — 1: 



(15) 



(16) 



A — 2A cos 



arg (2 + A 3 - 2VA 3 + l) 



(17) 



In figure [U we show the graphics of the positive solu- 
tion as a function of the parameter A. As expected, when 
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Figure 3: Plot of the general solution cr'/a for the cubic equa- 
tion l114[| as a function of the parameter A. 



A is negative (magnetization lying on the easy axis), the 
walls width is smaller than when A is positive (magneti- 
zation lying on the hard axis). We note that even when 
A > 0, i.e., the crystalline anisotropy does not favour 
the formation of walls, it is possible for the system to be 
divided in stable domains, as long as the dipolar inter- 
action is large enough compared to both the other two 
energy scales (as they appear in {Hp . 

An interesting situation that deserves a deeper anal- 
ysis is when < A < 1. Such case would occur, for 
example, if the energy scales satisfied J ~ g 2 fi%/a 3 and 
g 2 [i 2 B /a 3 3> AKa 3 . Then, we can expand (fT5|) and ob- 
tain: 



(^) = 2 2 / 3 + A + (A 2 ) , 



/ 2N + it{N-l) 



1/3 



in which there is an implicit dependence of cr'/a upon 
a. Although this procedure does not give exactly the 
same result as a numerical method, it prevents one from 
making further errors due to the integral that appears in 
e«(|,JV,p). 

Finally, it is interesting to study this model in the limit 
of d,D — > oo, where one expects to recover the result 
known as Landau-Lifshitz wall (see, for instance [23]). 
From equation (fTTj) , it is clear that the terms referring to 
the exchange and to the crystalline anisotropy are pro- 
portional to d~ l . The dipolar term has two parts: e z , 
which is due to the surface charge on the upper and lower 
faces of the slab, goes to zero as D — > oo; e x is just a fi- 
nite number for a constant aspect ratio p. Then, making 
p constant as d, D — > oo, it is clear that this term is 
proportional to d~~ 2 and is negligible if compared to the 
others. Hence, as pointed in [23j, the dipolar energy van- 
ishes in the limit of an "infinite" crystal and domain walls 
will be formed only if AK < 0. Minimizing the energy 
with respect to a and putting N = 2, we find the result: 



(!) 



J 



„> im \\2\AK\a 3 



which is just 6% greater than the value of the Landau- 
Lifshitz wall. 



from which we get, to first order in A: 



a 3 J 

9 2~ 



fAKa 6 



1/3 



/ 2N + tt(N~1) 
\l6y^^)e x (0,N,p) / 



V 9*1*% J ^e x (0,N,p) 



To avoid numerical problems and reduce the errors 
concerning the approximation lfl2|). it is useful to trans- 
form the general solution for the walls width in a self- 
consistent equation: 



(!) 



a 3 J 



1/3 



1 / 3 III. THEORETICAL HYSTERESIS LOOPS 



We now want to study how this system responds to an 

(18) external magnetic field applied along the x axis. It is ex- 
pected that, for strong enough fields, the magnetization 
will lie along the x axis and will follow the field direction. 
Our main objective is to determine the qualitative fea- 
tures of the hysteresis curves that would be observed for 
different slab widths. Therefore, we must generalize our 
previous model to include rotations of the domains mag- 
netization. Introducing the angle <j> between the mag- 
netization and the z axis, we can write its components 

(19) approximately as: 



M z (x) 



Mq cos . 



i=l 



i-iy 



erfc 



x-d{i- 1)/N 



V2a 



erfc 



di/N 



V2o 



6 



M x (x) = M„(l 



sin ( 



JV-l 



+ M 



sm< 



erfc 



V2a 



erfc 



x — d 
V2a 



(20) 



r 



assuming that the condition |(3]) is satisfied. In what 
follows, we will consider that the width of the walls is 
constant, to simplify the calculations. Using the formal- 
ism of last section, it is straightforward to conclude that 
an external magnetic field along the x axis usually will 
not substantially change the value of the minimum wall 
width. 

It is important to notice that this expression for M x , 
(|20"1) . as well as the previous one, |(4]), assumes that the 
domains are along the positive x direction. Although this 
feature does not influence any of the results obtained in 
the last section, it must be taken into account in this 
one, since we are dealing with spin rotation. If we were 
to make a complete description of this phenomenon, it 
would be necessary to include domains in the y direction, 
otherwise the magnetization inside the walls could never 
rotate properly. Hence, bidimensional domains would 
have to be considered, but this is beyond the scope of this 
article. Therefore, as the domain walls are very small, we 



will not treat their rotation, but only the domains rota- 
tion. This procedure will then be enough to give us the 
main qualitative characteristics of the hysteresis loops. 

Repeating the procedure of the last section, it is 
straightforward to calculate the total energy density re- 
ferring to the configuration ([20]), We substitute it in 
the expressions for the exchange, dipolar and crystalline 
anisotropic energies, obtaining the total energy accord- 
ing to (fTTjl . Leaving only the terms proportional to <j>, 
we obtain: 



E 
2M§ 



A (N, d, D) sin 2 



(21) 



B(N, d ,D) + §-c(N,^ 



sin < 



where H is the external magnetic field along the x 
direction and: 



J 



A (N, d, D) = 2np 



C 



-u 



« 

2 , 



U _u 

— he p 
P 



1 



2 sin u/2 /a\ /2sinu/2 

7T U 2 \d 



~4p I i 
'o 

J 



- 2 (f) 



7T U 
2 p-«/ 2 P 



sin 2 u/2 
sin 2 u/2N 



2 cos 



u(N-l) 
2N 



sin 7t 



2N 



COt (^) tan l2 1 ' 



du 



if 



sinh Qi) tan 2 (JL) [l - (-1)" cos u] du 



B (AT, d, D) = 2np 



e-" 2 (§) 



-2N + k{N - 1) + 2 
8y/n(d/a)(a/a) 

1 



o u 
a\ 1 sin u/2 



U _u 

he p 

P 



\dJ 2tt 



+2 

1 

2 



J 



5 2 M|/a E 



cot (^) tan 

n(N-l) 
8y/w(d/a)((r/a). 



(!) 



1-(-)V2tt(JV-1) 



r 



AKa 3 
2 \g 2 V 2 B /a 



a 
d 

- 1 



sin 2 u/2 
sin 2 u/2N 



2 cos 



u(N-l) 
2N 



sin ■ 



2N 



du 



(22) 



If the condition J3]) is satisfied, one can usually make 
the approximation C ~ 1/2, as long as there are few 
domains inside the slab. 

As we did not consider the terms independent of </) in 
the total energy expression lj2Tj) , we note that the energy 
of the configuration studied in the previous section, cor- 



responding to (f> = 0, would be E = 0. Indeed, the terms 
that appear in f2"Tj) are a combination of the exchange, 
dipolar and crystalline anisotropic energies referring to 
the components of the magnetization l(20|) that are not 
parallel to the z direction. It is clear that the total en- 
ergy will be non-zero only for configurations in which the 



7 



the magnetization is tilted by some angle with respect to 
the z axis. 

Now, it is possible to apply a procedure similar to the 
one developed by Stoner and Wohlfarth [24] to analyze 
the spin reorientation in the presence of a magnetic field. 
Minimizing equation lj22|) with respect to <j> leads us to 



dE 



B 



H 



Mn 



C cos < 



(23) 



2 A sin 4> cos <ft 

c)<p 

= 0, 

while the second derivative is given by 



|^ = 2Acos20 + fl! ■ -^-C : (24) 



positive magnetic field, the magnetization will lie along 
the positive x direction, since </>i = ir/2 will be the 
global minimum. As H diminishes, and eventually be- 
comes negative, condition i|26|) . at some moment, will no 
more be satisfied and two options will raise. If A < 0, 
while condition ([29]) will never be satisfied, condition lj27j) 
will be, and the system will jump to the configuration in 
which the magnetization lies along the negative x direc- 
tion (4> — — 7r/2). Hence, the upper curve of the hystere- 
sis loop has the shape of an abrupt step when A < 0. 
The remanent magnetization along the x axis is simply 



M r = M , 

and the coercive field is given by: 



(30) 



Equation (|23]l has two possible solutions: the first one 
refers to the two possible orientations for which the mag- 
netization lies on the x axis: 



cos^i = => 0i = ±- 



(25) 



and is a minimum as long as: 



if 4>\ = 7r/2 , or: 



(26) 



H, 



-Mn 



B-2A 
C 



(31) 



However, if A > 0, other stable states will raise before 
(f> = — ir/2 is reached, since condition |29|) is satisfied. 
Then, the system jumps continuously from minimum to 
minimum, as none of them can be a metastable state. 
This behaviour continues until H is strong enough to not 
satisfy condition |29]) anymore. Hence, condition (J27j) is 
satisfied and the magnetization lies along the negative 
x direction. In this case, in which A > 0, the upper 
curve has the shape of an inclined step. The remanent 
magnetization along the x axis is given by: 



24- - ( B + ^- Q C ! - 11 ■ 



if 4>i — —it/2. The second solution is: 



(27) 



Mr 



M 
Mq 



and the coercive field by 



if B > 2A 
otherwise 



(32) 



?2 = arcsm 



and it is a minimum for 



2A 



2A > 



B 



H 

Mo 



c 



2A 



(28) 



(29) 



Therefore, we can build a procedure to trace the upper 
curve of the hysteresis loop. Applying a strong enough 



"' = - U "«|?) 



(33) 



To obtain the lower curve of the hysteresis loop, we 
have to perform a slight modification. As previously ex- 
plained, the configuration considered assumes that the 
domain walls are along the positive x direction. However, 
after the magnetization is reoriented, we must consider 
a new configuration in which the domain walls are along 
the negative x direction, so we avoid additional errors for 
not considering domains along the y direction. Hence, we 
take: 



J 



M x {x) = Mo (-1 



N-l 

i=i 



-^71 



Mn 



erfc 



V2a 



erfc 



V2a 



(34) 
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(35) 



Equation <f2T|) is then replaced by 
I 



-2 = A (N, d, D) sin 2 + 



Following the same procedure as before, we obtain that 
the lower curve has the same shape of the upper one, 
and the remanent magnetization along the x axis and 
the coercive field are symmetrical to the ones presented 
previously, {30} - {33]) , 

Therefore, we conclude that, if A < 0, the hysteresis 
loop is squared, the remanent magnetization is Mq and 
the coercive field is given by lf3l]) . If A > 0, instead, the 
hysteresis loop is inclined, the remanent magnetization 
is given by 11321) and the coercive field by (|33|) , 



IV. APPLICATION TO MnAs : GaAs 

Several works were done regarding magnetic properties 
of MnAs thin films grown over GaAs substrates. Here we 
will discuss possible applications of our model to the ex- 
perimental hysteresis loops recently obtained [27l.l28l.l29l]. 
To make correspondence between our model and the real 
system, it is useful to identify the following crystallo- 
graphic directions of MnAs to the axis considered in fig- 
ure 1: [1120] = x, [0001] = y and [1100] = z. 

First of all, let us outline the parameters of the MnAs 
thin films that we will use in our calculations. As mea- 
sured in [HO], the anisotropy constants associated to the 
z and to the x axis are, respectively, K z — 7 ■ 10 6 erg/cm 3 
and K x = 7.4 • 10 6 erg/cm 3 . Hence, the x axis is the easy 
one, and we can use AK = 0.4 • 10 6 erg/cm 3 . From the 
crystalline structure of MnAs Q, which is in fact hexago- 
nal, we estimate the value of the equivalent cubic lattice 
parameter to be a — 4 A. The gyromagnetic factor was 
estimated as g = 4.5, from which we obtain a magneti- 
zation of Mq = 0.65 • 10 6 A/m, that is very close to the 
experimental value measured of Mo — 0.67- 10 6 A/m [3lj . 
As we are interested in orders of magnitude, we estimate 
J = 4.5 meV from the Curie temperature of MnAs. Fi- 
nally, as these films have constant thickness, we take the 
fixed value D — 130 nm, which is the thickness of the 
sample used in [29] and of the same order of magnitude 
of usual samples. 

Magnetic Force Microscopy (MFM) images of 
MnAs:GaAs films have suggested that [li, [33], as the 
temperature rises in the coexistence region (implying in 
smaller widths of the ferromagnetic stripe), the system 
undergoes a transition from a configuration in which the 
magnetization lies along the easy x axis to another one 
in which it lies along the growth axis (the hard z axis) 
and is divided in three domains. Let us verify the predic- 



r 



AE 




Figure 4: Plot of the energy difference between the configura- 
tion in which the magnetization lies along the x axis and the 
one in which there are three domains lying along the z axis as 
a function of the aspect ratio p of the ferromagnetic stripe. 



tions of our model for these situation: substituting the 
experimental parameters of MnAs in the self-consistent 
equation lfl9|) for N = 3 domains, we obtain that the 
domain walls are stable and that their width is given by 
a = 3.8a. The solution to the cubic equation l[T4"ll used 
was l|18p , since the experimental parameters of MnAs im- 
ply A< 1. A numerical calculation in fact gives a width 
of about 5a, what confirms our expectations that the ap- 
proximative self-consistent method gives errors less than 
20%. 

Using equation (J2TJ) , it is possible to obtain the en- 
ergy difference between the two configurations seen in 
the MFM images. In figure HJ we show the plot of this 
difference as a function of the aspect ratio p and notice 
that there is a transition at p c ~ 1.5. The measurements 
realized by Coelho et. al. [29j suggest that this transition 
takes place around p c « 2.9. However, as mentioned pre- 
viously, we do not take into account the domains along 
the y direction or the inter-stripe interaction, a feature 
always observed in the MFM images. This seems to be 
particular important not only to provide more precise 
values for p c but also to ensure that the most stable con- 
figuration has 3 domains. Without taking into account 
such features, our model would predict also transitions to 
N > 3 domain configurations before the one at N = 3, 
what is not observed in the MFM images. Hence, we can 
only predict a qualitative behaviour of the real system. 
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Figure 5: Plot of the hysteresis loop predicted by the model 
for p = 1.9. The vertical axis represents the relative magneti- 
zation M/Mq along the x axis and the horizontal axis repre- 
sents the external magnetic field H applied in the x direction 
in Oe. 
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Figure 6: Plot of the hysteresis loop predicted by the model 
for p — 1.7. The vertical axis represents the relative magneti- 
zation M/Mo along the x axis and the horizontal axis repre- 
sents the external magnetic field H applied in the x direction 
in Oe. 



It is interesting to plot the theoretical hysteresis loops 
predicted by our model for a N = 3 domain configuration 
with the MnAs experimental parameters. Applying the 
procedure developed in the previous section, we see that 
B > for any value of the aspect ratio p and A is negative 
until p « 1.8 , where it becomes positive. Hence, it is 
expected that the hysteresis shape is squared until p w 
1.8, where it becomes inclined. For 1.65 < p < 1.8, the 
remanent magnetization is expected to be yet Mo, and 
only for p < 1.65 it will start to decrease. Figures [5j 
[6] and [Jj show the three different possible shapes of the 
hysteresis loops, according to the procedure of the last 
section. 

It is worth to notice that hysteresis loops of shapes 
similar to figures [71 and [6] (but more rounded) were ob- 
served in [28| and [29j , respectively, while shapes like the 
one in figure 4 were seen in both of them. Coelho et. 
al. showed that the change in the shape of the hysteresis 
loops occurs for p ~ 2.9, which is not too far from our 
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Figure 7: Plot of the hysteresis loop predicted by the model 
for p = 1.5. The vertical axis represents the relative magneti- 
zation M/Mo along the x axis and the horizontal axis repre- 
sents the external magnetic field H applied in the x direction 
in Oe. 
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Figure 8: Plot of the (a) modulus H c of the coercive field (in 
Oe) and (b) of the relative remanent magnetization M/Mo 
along the x axis as a function of the aspect ratio p of the 
ferromagnetic stripe. 



prediction (p « 1.8). However, our coercive fields are 
one order of magnitude greater than the values obtained 
by both [2§| and (28|. It is interesting to notice that 
simulations performed by Engel-Herbert et. al. jH| con- 
sidering other configurations and other methods led to 
the same order of magnitude for the coercive fields than 
our model. The graphics of the coercive field and the re- 
manent magnetization predicted by our model are shown 
in figure where it is easy to verify the paths for the 
transitions among the three different types of hysteresis 
loops pointed before. 

Finally, we point out that, in [23) Takagaki et. al. 
found a hysteresis loop similar to the squared one (figure 
[5]) for MnAs thin films and another similar to the in- 
clined one (figure [6]) (but more rounded, again) for disks 
of MnAs fabricated from thin films. It is evident that our 
model cannot be applied to disks (in such case, there is 
the possibility of more complex configurations, like vor- 
tices, for example), but it may give a hint about their 
physical behaviour. 
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V. CONCLUSIONS 

In this article, we have shown, using an approxima- 
tive analytical self-consistent equation, that multiple- 
domains configurations with sharp walls in a ferromag- 
netic slab can be stable. We arrived at ratios among 
the three typical magnetic energy scales of the system 
(exchange, dipolar and crystalline anisotropic) that can 
determine the stability of these configurations and also 
give the order of magnitude for the walls width. Although 
the approximations done to achieve these equations intro- 
duces imprecision on the predictions of the model, they 
do not significantly change the orders of magnitude in- 
volved. Moreover, this model shows that, even when the 
crystalline anisotropy prefers spread walls, the dipolar 
interaction can compensate it to form sharp ones. In 
what concerns thin films of MnAs grown over GaAs sub- 
strates, which could be an observable realization of our 
model, we corroborated the suggestions based on MFM 
images that predicts the formation of three-domain con- 
figurations along the hard axis, for temperatures above 
25 °C. The transition between this state and the config- 
uration in which the magnetization lies completely along 
the easy axis was also predicted, but the value of the film 
aspect ratio for which this transition occurs was far from 
the experimental one. 

In addition, we compared the hysteresis loops that ap- 
pears when an external magnetic field is applied along the 
easy axis direction for such three-domain configurations. 
Qualitatively, we obtained, using an approach similar to 
Stoner and Wohlfarth , all the three shapes of loops 
observed in the literature (squared, inclined with large re- 
manent magnetization and inclined with small remanent 



magnetization). The main differences are that the in- 
clined loops measured experimentally are more rounded 
than ours and that the experimental coercive fields are 
one order of magnitude smaller. 

Several factors help us to understand why all these 
differences between the predictions of the model and the 
experimental measurements occur. Firstly, the simple 
model we used here does not consider features that are 
essential in the MnAs:GaAs real system, like the modu- 
lation along the y direction and the inter-stripe dipolar 
interaction. To include these features, a more sophisti- 
cated model with bidimensional domains and topological 
defects would be necessary. Moreover, in what concerns 
the hysteresis loops, the Stoner- Wohlfarth method deals 
only with the collective spin rotation, and does not take 
into account nucleation or pinning, which can be respon- 
sible for the rounded shape of the curves and the smaller 
coercive fields observed in the experiments. As already 
pointed out in due to the lack of complete knowledge 
about the exact geometrical forms of the ferromagnetic 
stripes, the materials inhomogeneities (that can induce 
nucleation) and the correct microscopic parameters, we 
cannot expect an exact reproduction of the experimen- 
tal hysteresis loops. Nonetheless the qualitative physical 
properties predicted here give an insight for a more com- 
plete understanding of the complex domain structure of 
systems like MnAs:GaAs films. 
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